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Moist  air  is  characterized  Tor  the  frequency  range  1-1000  GHz  as  a  nonturbulent  propagation 
medium  described  by  meteorological  parameters.  An  adequate  spectroscopic  data  base  for  air  consists 
of  three  terms:  (1)  resonance  information  for  29  H20  lines  up  to  1097  GHz  and  44  02  lines  up 
to  834  GHz  in  the  form  of  intensity  coefficients  and  center  frequency  for  each  line;  (2)  an  empirical 
water  vapor  continuum  spectrum;  and  (3)  a  liquid  water  attenuation  term  for  haze  and  cloud  conditions. 
This  data  base  is  the  heart  of  two  computer  programs  which  calculate  and  plot  attenuation  rates 
(in  decibels  per  kilometer),  refractivity  (in  parts  per  million),  and  refractive  dispersion  (in  parts 
per  million).  The  first  covers  the  troposphere  and  requires  pressure,  temperature,  and  relative 
humidity  as  input  data.  The  second  addresses  isolated  line  behavior  in  the  mesosphere  wherein 
the  geomagnetic  field  strength  H  is  an  additional  input  parameter  due  to  the  Zeeman  effect  of 
the  02  molecules.  Each  oxygen  line  splits  proportionally  with  H  into  numerous  sublines,  which 
are  juxtaposed  to  form  Zeeman  patterns  spread  over  a  megahertz  scale.  Patterns  of  three  main 
polarization  cases  are  considered.  Various  typical  examples  for  a  model  atmosphere  demonstrate 
the  utility  of  the  approach,  provide  new  information,  and  underline  the  serious  role  that  water 
vapor  plays  above  120  GHz. 


1.  INTRODUCTION 

The  physical  properties  of  the  neutral  atmosphere 
influence  radio  wave  propagation.  For  the  fre¬ 
quency  range  1-1000  GHz,  inherent  losses  are  due 
to  spectral  absorption  by  the  molecules  H20  and 
02.  Further  propagation  limitations  are  caused  by 
turbulence  and  atmospheric  layering  effects,  and 
>hy  thermal  noise  emitted  from  the  absorbing  air 
Qmass.  Intermittently,  of  course,  rain  is  the  dominant 
CI3bsorber,  at  least  in  the  15-  to  500-GHz  range 
Cjtrane,  1981], 

With  the  current  high  interest  in  millimeter  and 
"■sjibmillimeter  waves,  there  is  need  for  a  reliable 
— model  to  predict  average  loss  and  delay  effects 
^"‘Ttom  easy-to-obtain  meteorological  data.  Such  a 
—  model  would  find  considerable  practical  application 
■Smrough  conversion  of  basic  climatological  variables 
EZE.e.,  barometric  pressure  P,  temperature  T,  relative 
“timidity  RH)  into  transfer  characteristics  of  a  radio 
path.  The  frequency  of  1000  GHz  is  a  kind  of  new 
frontier  for  a  growing  technology  encompassing 
microwave  and  optical  principles.  Millimeter  wave 
systems  are  attractive  because  of  their  ability  to 
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penetrate  a  somewhat  opaque  atmosphere  (haze, 
fog,  clouds,  dust,  smoke,  light  rain)  under  circum¬ 
stances  in  which  electro-optical  and  infrared  sys¬ 
tems  normally  fail.  Accurate  and  detailed  knowledge 
of  atmospheric  transmission  properties  is  essential 
in  evaluating  the  advantages  that  millimeter  waves 
might  have  over  shorter  wavelengths. 

Water  in  both  vapor  and  liquid  states  is  the  major 
deterrent  to  an  unrestricted  exploitation  of  millime¬ 
ter  and,  more  so,  of  infrared  wavelengths.  For  most 
applications,  the  operation  of  ground-based  systems 
is  limited  in  frequency  to  window  regions  W\-W5, 
these  being  the  gaps  between  molecular  absorption 
lines  and  bands  as  schematically  given  in  Table  1. 

Mean  values  and  limits  of  the  following  propaga¬ 
tion  effects:  absorptive  loss  of  coherent  radiation, 
time  of  propagation  between  two  points  (both  of 
which  are  discussed  in  this  paper),  refractive  ray 
bending,  generation  of  incoherent  noise,  and  scin-  ' 
filiations  due  to  random  fluctuations  of  the  medium 
in  space  and  time,  can  be  considered  for  modeling. 
The  characteristics  of  a  short,  horizontal  radio  path 
are  approximately  modeled  by  one  value  of  an 
average  complex  refractivity  N.  Cumulative  behav¬ 
ior  (e.g.,  a  ground- to-satellite  path)  is  calculated 
by  assuming  a  spherically  stratified  atmosphere  in 
which  each  layer  has  a  constant  N  value.  This 
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TABLE  I .  Millimeter  wave  window  ranges  of  the  atmosphere. 


Absorption  feature 

Window  range 
(GHz) 

Attenuation  range 
at  sea  level 
(dB/km) 

22-GHz  H.O  line 

W\ 

24-48 

0. 1-0.3 

60-GHz  O,  line  complex 
W2 

70-110 

0.3-1 

1 19-GHz  O,  line 

m 

120-155 

1-2.5 

183-GHz  H20  line 

WA 

190-300 

2-10 

325-GHz  H,0  line 

W5 

335-355 

5-20 

380-GHz  H.O  line 
(and  1823  more  lines 
of  the  rotational 
band  up  to  25  THz) 

layer-by-layer  method  uses  numerical  integration 
techniques.  Refinements  consider  refractive  and 
dispersive  ray  bending  at  low  elevation  angles 
[Hopponen,  1980] . 

Since  the  pioneering  work  by  Van  Vleck  [1947] 
and  Birnbaum  [  1953]  much  progress  has  been  made 
in  describing  atmospheric  millimeter  wave  spectra. 
Recent  advances  are  in  the  spectroscopic  data  base 
[Liebe  et  al.,  1977;  Poynter  and  Pickett,  1980; 
Rothman,  1981] ,  in  line  shape  theory  [ Rosenkranz , 
1975;  Lam,  1977;  Smith  and  Guiraud,  1979;  Smith, 
1981] ,  in  measurements  of  refraction  [ Kemp  et  al., 
1978] ,  in  modeling  [Falcone  et  al.,  1979;  Hopponen, 
1980;  Hill  and  Clifford,  1981] ,  in  reposing  anoma¬ 
lous  window  absorption  [Emery  et  al.,  1980],  in 
collecting  available  field  data  to  support  an  empirical 
water  vapor  continuum  [Crane,  1980] ,  in  reviewing 
window  excess  absorption  [Liebe,  1980] ,  in  treating 
liquid  water  uptake  in  clear-to-hazy  air  at  high 
relative  humidities  [Nilsson,  1979],  and  in  the 
submillimeter  wave  dielectric  constant  for  water 
[Simpson  et  al.,  1979].  A  unified  model  of  these 
new  results  is  given  in  two  computer  programs; 
Program  PI  covers  attenuation  rates  a(/)  and  dis¬ 
persion  D(f)  in  the  height  range  h  =  0-30  km; 
Program  P2  calculates  isolated  line  behavior  espe¬ 
cially  Zeeman  patterns  of  oxygen  lines  in  the  meso- 
phere  ( h  =  30-100  km).  Both  programs  are  formu¬ 
lated  in  engineering  terms  and  are  fully  traceable, 
but  void  of  quantum-mechanical  complications. 
They  can  serve  as  a  reference  for  comparisons  with 
future  data,  which  are  presently  very  scarce  for 
the  200-  to  1000-GHz  range.  In  this  range,  the 


computation  rests  foremost  on  theory,  while  the 
confidence  limit  below  100  GHz  is  estimated  to 
be  better  than  4%,  based  on  extensive  laboratory 
testing. 

In  the  past,  computer  programs  were  written 
[Falcone  et  al.,  1979;  Harries,  1980;  Hill  et  al., 
1980]  with  spectral  line  data  drawn  from  the  Air 
Force  Geophysics  Laboratory  tape  ( Rothman 
[1981]  and  Rothman  et  al.  [1981]  ;  these  compila¬ 
tions  are  available  as  the  main  atlas  of  some  159,000 
transitions  and  a  minor  constituents  atlas  in  the 
form  of  magnetic  tapes  from  the  National  Climatic 
Center  of  the  National  Oceanic  and  Atmospheric 
Administration,  Digital  Product  Section,  Federal 
Building,  Asheville,  North  Carolina).  The  sheer 
number  of  lines  (>105)  makes  such  calculations 
quite  time-consuming,  nontransparent,  and  costly. 
Up  to  1000  GHz,  the  two  principal  absorbing  species 
are  H20  and  02,  and  simplifications  can  be  em¬ 
ployed  such  as  (1)  an  intensity  cutoff  (>2  X  10~5 
dB/km  at  line  center)  to  reduce  the  number  of 
contributing  lines  and  (2)  an  adjustable  resolution 
factor  determining  the  frequency  grid  for  plotting, 
which  still  captures  all  extremes  without  an  exces¬ 
sive  number  of  steps.  On  the  other  hand,  important 
additions  are  the  overlap  correction  for  the  60-GHz 
02  band  and  the  Zeeman  splitting  matrices  for 
mesospheric  02  lines. 

2.  REFRACTIVITY  OF  MOIST  AIR 

Amplitude  and  phase  response  of  a  planar  radio 
wave  traveling  the  distance  L  and  starting  with  the 
field  strength  E0  is  described  by 

E  =  E0  exp  (TL)  =  E0exp  [,/0.02096/(10‘  +  N)L\  (I) 

Frequency  f  is  in  gigahertz.  (GHz)  throughout  the 
paper  and  T  is  the  propagation  constant. 

Germane  to  any  propagation  mode)  is  a  conve¬ 
nient  macroscopic  measure  of  the  interaction  be¬ 
tween  radiation  and  the  absorbing  species  in  moist 
air.  Complex  refractivity  N  (in  parts  per  million), 
expressed  in  terms  of  measurable  quantities,  pro¬ 
vides  that  role.  For  air,  N  consists  of  three  compo¬ 
nents 

N  =  N„+  D(f)  +  jN”(f)  ppm  (2) 

namely,  the  frequency  independent  refractivity  N0 
plus  various  spectra  of  refractive  dispersion  D(f) 
and  absorption  N"(f).  Usually,  the  imaginary  part 
of  (2)  is  expressed  as  the  specific  power  attenuation 


ATTENUATION  AND  PHASE  OF  RADIO  WAVES  1185 

a  and  the  real  part  determines  the  phase  delay  p 
(with  reference  to  vacuum);  that  is, 

a  =  0.  1820/jV"  dB/km 

(3a) 

3  =  0.02096/(lVo  +  D)  radians/km 

Accordingly,  the  propagation  constant  r  and  the 
excess  propagation  delay  time  t  are 

T  = -0.1151a +7(2.096  x  104/ +  3)  1/km 

,  (3*) 

t  =  (3/2it/)103  =  3.336(W0  +  D)  ps/km 

An  overview  of  notation  used  to  formulate  a(N ") 
and  3(1V0,  D)  is  given  in  Table  2. 

The  physical  state  of  moist  air,  which  determines 
the  absorber  population  in  a  radio  path,  is  described 
by 

p  dry  air  pressure  (barometric  pressure  P  =  p  +  e), 
kPa  (1  kPa  =  10  mbar); 

0  relative  inverse  temperature,  equal  to  300/  T  (T  in 
°K.); 

e  water  vapor  partial  pressure,  kPa; 
tv  liquid  water  concentration,  g/m3. 

The  number  M  of  O  (20.845%  in  dry  air)  and 
H20  molecules  per  unit  air  volume  is  given  by  the 
ideal  gas  law: 


TABLE  2.  Notation  overview. 


Notation 

Unit  of  measurement 

Propagation  parameters 

Frequency 

/(I) 

GHz 

Attenuations 

a  (3),  a.  (22).  a,  (25).  a,  ,  (28).  a,  (29) 

dB/km 

Refractivities,  dispersion 

N  (1).  N„  (10).  N"  (11),  AT  (20),  (21).  D  (12) 

ppm 

Orientation  angle 

0(29) 

Attenuation  response 

5  (Figure  9) 

% 

Physical  measurables 

Height 

A  (9) 

km 

Pressure,  temperature,  relative  humidity 

P(h).  T(h).  RH  (h)(9) 

kPa.  K..  % 

Magnetic  field  strength 

H  (26) 

G 

Physical  variables 

Molecular  number  density,  relative  mass 

M  (4),  m  (23) 

cm  ’.  1 

Relative  temperature,  dew  point 

«  (4).  a,  (5) 

1 

Partial  pressures 

p,  e  (4).  e 

kPa 

Droplet  concentration 

w<7) 

g/m3 

Spectroscopic  parameters 

Line  shapes,  strength 

F"  (13).  r  (14).  S  (17) 

GHz ■*.  MHz 

Center  frequency,  widths,  overlap 

v„  (13),  y  (17),  y„  (23).  0(17) 

GHz.  GHz.  1 

Line  coefficients 

*,  ,<17) 

(see  Table  4) 

Zeeman  shift  and  relative  strength 

ri  (26).  4  (27) 

1 

Transition  labels 

ir,  o'  (Table  6) 

O,  quantum  numbers 

K.  M  (Table  6) 

Complex  dielectric  constant  of  water 

r  =  f  ’  +  jr"  (22 a) 

1 

Number  in  parentheses  denotes  equation  where  symbol  is  defined. 


M(Oj  +  HjO)  =  (0.5034 p  +  2.415e)d  1017  cm"3  (4) 

Thus  p-e-0  combinations  account  for  the  number 
of  molecular  absorbers  in  moist  air.  In  the  water 
vapor  state,  the  saturation  pressure  es  at  a  given 
‘temperature’  0,  (dew  point)  cannot  be  exceeded; 
that  is  [Bogel,  1977], 

e  £  es  =  2.40903  X  j0|lo-,  8548i>  kPa  (5) 

e,  “  3.5018 

(Equation  (5)  describes  the  data  of  the  standard 
Smithsonian  Meteorological  Tables  to  within  ±0.2% 
for  the  range  ±40°C.)  The  vapor  pressure  with 
respect  to  saturation  defines  relative  humidity  as 

RH  =  (e/eJlOOg  100%  (6) 

The  humidity  condition  RH  =  100%  is  a  delicate 
balance  point  for  physical  phase  changes.  A  drop 
of  1  K  decreases  the  water  vapor  in  saturated  air 
between  5%  (40°C)  and  8%  (— 10°C).  The  excess 
water  vapor  pressure  ex  =  e(0,)  -  es( 0*  >  0,) 
is  converted  into  a  liquid  water  droplet  concentra¬ 
tion 

w  =  7.219etdf  g/m3  (7) 

The  conversion  process,  however,  starts  well  below 
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TABLE  3.  The  mc.i-1  atmosphere  (U  S.  Standard  Atmosphere,  1976). 


Height,  h 
(km) 

Dry  air 
pressure,  p 

Temperature.  T 

<°K) 

Maxi¬ 

mum 

relative 

humidity 

(/?«)„ 

(%) 

Magnetic 

field 

strength  H 
(G) 

Program 

Temperature 

parameter, 

0  =  300/ T 

Saturation 
pressure,  e, 
(equation  (5)) 
(Pa) 

0 

101.325  kPa 

288.15 

100 

PI 

1.041 

1706 

1.6 

83,527  It  Pa 

277.75 

100 

1.080 

847.7 

3 

70.121  kPa 

268.66 

100 

1.117 

434.0 

6 

47.217  kPa 

249. 19 

100* 

1-1000  GHz 

1.204 

88.1 

10 

26.499  kPa 

223.25 

60 

1.344 

6.41 

16 

10.352  kPa 

216.65 

3 

1.385 

(172)+ 

20 

5.529  kPa 

216.65 

2 

1.385 

<1  72)+ 

30 

1.197  kPa 

226.51 

1 

1.324 

9.36 

30 

1197.000  Pa 

226.51 

0.2-0. 8 

P2 

1.324 

40 

287.140  Pa 

250.35 

0.2-0. 8 

1. 198 

50 

79.779  Pa 

270.65 

0.2-0. 8 

1.108 

60 

21.95 8  Pa 

247.02 

0.2-0. 8 

1.2)5 

70 

5.221  Pa 

219.59 

0.2-0. 8 

MO,) 

1.366 

80 

1.052  Pa 

198.64 

0.2-0. 8 

±50  MHz 

1.510 

90 

0.184  Pa 

187.21 

0.2-0. 8 

1.603 

100 

0.032  Pa 

198.99 

0.2-0.8 

1.508 

•Height  profiles  of  humidity  above  h  =  6  km  are  discussed  by  Ellsaesser  el  al.  ( 1980] . 
tValue  over  ice. 


RH  =  100%  due  to  water  uptake  by  suspended, 
invisible  aerosol  particles  [Nilsson,  1979],  The 
liquid  water  concentration  in  ‘clear’  air  with  RH 
>  70%  can  range  somewhere  between 

w(RH)  =  0.001  and  0.1  g/m1  (8) 

The  objective  at  this  point  is  to  express  the 
quantities  a(/),  D(f),  and  N„  in  terms  of  p-e-ft-w 
data,  subject  to  the  following  conditions: 

Frequency 


N„  =  2.589pS  +  41.6e02  +  2.39eO  ppm  (10) 

and  can  be  calculated  in  a  straightforward  manner. 
Water  vapor  refractivity  is  about  16  times  more 
effective,  on  a  per  molecule  basis,  than  dry  air 
in  generating  propagation  phenomena  such  as  delay, 
ray  bending,  ducting,  scintillations,  etc. 

The  absorption  and  dispersion  spectra  are  formu¬ 
lated  from  line  contributions,  a  continuum  N"c,  and 
a  liquid  water  extinction  AT";  i.e., 

N"(f)  =  £  (SF" ),  +  N"  +  N".  ppm  (11) 


/=  1-1000  GHz 


Altitude 


and 


h  =  0-t00km 

Relative  humidity 

RH  =  0-100% 


(9) 


Magnetic  field  strength 

H  =  0. 1-0.8  gauss  (G) 

Numerical  and  graphical  examples  are  then  calcu¬ 
lated  for  the  model  atmosphere  specified  in  Table 
3,  using  various  combinations  of  (9). 

The  frequency-independent  refractivity  is  given 
by  [Liebe  el  al.,  1977] 


D{f)=^(SF'),  ppm  (12) 

Both  expressions  require  more  elaboration,  as  de¬ 
tailed  below. 

First,  the  line  spectra  of  absorption  SF ”  and  of 
refractive  dispersion  SF'  have  strength  S  in  units 
of  kilohertz  and  shape  factors  F'  and  F"  in  units 
of  (GHz)  '.  The  summations  consider  j  =  44  02 
plus  29  H20  lines  (line  spectra  of  the  trace  gases 
O,,  CO,  N20,  SOj,  NH}.  etc.  [Waters,  1976; 
Poynter  and  Pickett,  1980;  Rothman,  1981;  Rothman 
et  al.,  1981]  are  neglected).  Common  to  each  line 


is  a  pair  of  intensity-versus-frequency  distribution 
functions  [e.g.,  Waters,  1976;  Kemp  et  al.,  1978]. 
In  the  atmosphere,  the  shape  factors  named  after 
Van  Vleck-Weisskopf  (VVW)_and  Gross  are  used. 
We  chose  the  VVW  shape  as  modified  by  Ro- 
senkranz  [1975]: 


,  //\r7-<»o 

\l’o/  -(>'(.  -/) 


~/)d  7  -  (k0  +f)d 

+  y  +K  +/)2  +  y2 


(13) 


and 


F,  _  .  K  -/)  +  yd  +  (f0  +/)  +  yd 


(v0-f)2  +  y2  (‘'o  +f)2  +  y 


-3 


(14) 


where  v0  is  the  molecular  line  center  frequency, 
7  is  the  width,  and  d  <  ±  1  is  the  overlap  interference. 
Above  heights  of  about  15  km,  isolated  line  behavior 
takes  over,  limiting  each  molecular  resonance  to 
a  megahertz  scale  frequency  span.  For  isolated  line 
calculations,  (13)  and  (14)  can  be  replaced  by 
Lorentzian  shapes 


K  =  y/  [(vo  -f)2  +  y2] 

F'l  =  (v„  -/)/  [(V„  -f)2  +  73] 


(15) 


Maximum  absorption  at  /  =  v„  and  peak  dispersion 
at  /  =  r »  +  y  are  given  by 


N"m  =  S/7  ±D„  =  S/2-y  ppm 


(16) 


The  line  parameters  are  calculated  as  follows 
[ Liebe  and  Gimmestad,  1978] : 


02  in  air 


S,  kHz 
7,  GHz 


atp9’  exp  ]a2(l  -  0)1 
a,(p  +  l.3*)6°‘' 
aiPV' 


(17) 


H20  in  air 


S,  kHz 
7,  GHz 

a 


b,e§  exp  [62(1  - 
6,(4.80e  +  p)«°‘ 
0 


An  invariable  line-data  base  (v0  and  the  spectro¬ 
scopic  coeff  dents  a,  , ,  613 )  is  used  in  conj unction 
with  meteorological  conditions  expressed  through 
sets  of  p-d-e  data.  The  line  base  is  derived  from 
experiments  [Liebe  et  al.,  1977;  Kemp  et  al.,  1978] 
and  spectroscopic  parameter  compilations  [Poynter 
and  Pickett,  1980;  Rothman,  1981]  and  is  listed 
in  Table  4.  Line  intensities  (equation  (16))  of  a  few 
selected  lines  are  independent  of  temperature  when 
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T=  273 

150  K 

6  =  1.1 

2.0 

a2  =  2.1  In  6/(6  - 

1) 

2.00 

1.46 

62  =  2.9  In  6/(6  - 

1) 

2.76 

2.02 

(18) 

Lines  with  a2,  b2  coefficients  lower  than  those 
prescribed  by  (18)  increase  in  intensity  when  the 
temperature  drops  and  vice  versa. 

Dispersion  D(f)  addresses  only  frequency  de¬ 
pendent  contributions.  The  term  2/v0  is  the  finite 
value  of  the  first  two  terms  in  the  shape  factor 
F’(f)  when /  — »  0  and,  consequently,  is  subtracted 
from  (14).  Without  this  measure,  for  example,  the 
sum  of  all  rotational  H20  lines  (that  is,  1838  lines 
between  22  GHz  and  31  THz)  is  close  (—2.5%) 
to  the  measured  dipole  orientation  term  41.6eb2 
in  (10)  [Hill  and  Clifford,  1981],  while  the  lines 
listed  in  Table  4  yield  a  residual  contribution  fitted 
by 

=  X  [SF'(/— .0)], 

=  (4.75e  +  1.14  x  10  >)«5  ppm  (19) 

Second,  the  continuum  spectrum  N"c  ~  N"  + 
N"d  consists  of  N",  the  far- wing  absorption  by  very 
strong  infrared  H20  lines  (plus  additional  contribu¬ 
tions  which  are  not  fully  understood)  and  N"d,  the 
nonresonant  dry  air  spectrum.  We  opted  for  the 
empirical  expression  [  Waters,  1976] 

N”  =  1.9/ep-d’1  x  10  ‘  ppm  (20a) 

which  fitted  available  field  data  reasonably  well 
[Crane,  1980], 

The  magnitude  of  H”  is  about  a  factor  of  5  above 
the  far-wing  contribution  (first  term  of  (206)),  which 
is  computed  when  the  complete  rotational  H20 
spectrum  is  taken  into  account.  The  excess  absorp¬ 
tion,  so  far,  has  not  found  a  generally  acceptable 
explanation.  An  earlier  laboratory  experiment  by 
us  yielded  (for  f  <  100  GHz) 

N",  =  (0.24/l5e/>«”  +  2.2/' V«IJ)  x  10  *  ppm 

(206) 

The  correct  magnitude  of  N"  is  a  major  factor  in 
the  window  ranges  W2-W5.  In  view  of  the  conflict¬ 
ing  evidence  [Liebe,  1980] ,  the  alternatives  of  either 
(20a)  or  (206)  are  not  too  different  (see  Figure  2 
in  section  3).  Various  researchers  have  conjectured 
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TABLE  4a.  Date  base  for  02  (ara,)  tad  H30  (b,-b,  )  spectral  lines  in  air. 

Center 
frequency.  v0 

Strength,  a, 

Temperature 

Width,  a, 

Interference,  a. 

Temperature 

Quantum  number 

(GHz) 

(ItHz/kP a) 

exponent,  a2 

(GHz/kPa) 

(1/kPa) 

exponent,  as 

identification,  K‘ 

0" 

0.307  E-3  (ppm/ltPa) 

0.12  (0.056)  E- 

1 

1*  to  37* 

50.47381 

0.94  E-6 

0.969  E+l 

0.86  E-2 

0.520  E-2 

0.179  E+l 

37 

50.98742 

0.244  £-5 

0.869  E+l 

0.87  E-2 

0.550  E-2 

0.169  E+l 

35 

51.50311 

0.604  E-5 

0.774  E  +  l 

0.89  E-2 

0.560  E-2 

0.177  E+l 

33 

52.02124 

0.141  E-4 

0.684  E+l 

0.92  E-2 

0.550  E-2 

0.181  E  +  l 

31 

52.54227 

0.308  E-4 

0.600  E+l 

0.94  E-2 

0.569  E-2 

0.179  E+l 

29 

53.06683 

0.637  E-4 

0.522  E  +  l 

0.97  E-2 

0.528  E-2 

0.189  E+l 

27 

53.59570 

0.124  E-3 

0.448  E  +  l 

0.100  E-l 

0.544  E-2 

0.183  E+l 

25 

54.12997 

0.2265  E-3 

0.381  E  +  l 

0.102  E-l 

0.480  E-2 

0.199  E+l 

23 

54.67115 

0.3893  E-3 

0.319  E+l 

0.105  E-l 

0.484  E-2 

0.190  E+l 

21 

55.22137 

0.6274  E-3 

0.262  E+l 

0.1079  E-l 

0.4)7  E-2 

0.207  E+l 

19 

55.78382 

0.9471  E-3 

0.212  E+l 

0.1110  E-l 

0.375  E-2 

0.207  E+l 

17 

56.26477 

0.5453  E-3 

0.100  E-l 

0.1646  E-l 

0.774  E-2 

0.890 

D'*  ''  ,5 

56.36339 

0.1335  E-2 

0.166  E+l 

0.1144  E-l 

0.297  E-2 

0.229  E+l 

56.96818 

0.1752  E-2 

0.126  E+l 

0.1181  E-l 

0.212  E-2 

0.253  E+l 

13 

57.61249 

0.2125  E-2 

0.910 

0.1221  E-l 

0.940  E-3 

0.376  E+l 

11 

58.32389 

0.2369  E-2 

0.621 

0.1266  E-l 

-0.550  E-3 

-0.111  E+2 

D2  r  9 

58.44658 

0.1447  E-2 

0.827  E-l 

0.1449  E-l 

0.597  E-2 

0.790 

59.16422 

0.2387  E-2 

0.386 

0.1319  E-l 

-0  244  E-2 

0.700  E-l 

7 

59.59098 

0.2097  E-2 

0.207 

0.1360  E-l 

0.344  E-2 

0.490 

5' 

60.30604 

0.2109  E-2 

0.207 

0.1382  E-l 

-0.435  E-2' 

0.680 

D3  r  5 

60.43478 

0.2444  E-2 

0.386 

0.1297  E-l 

0.132  E-2 

-0.120  E+l 

61.15057 

0.2486  E-2 

0.621 

0.1248  E-l 

-0.360  E-3 

0.584  E+l 

9’ 

61.80016 

0.2281  E-2 

0.910 

0.1207  E-l 

-0.159  E-2 

0.286  E+l 

11* 

62.41122 

0.1919  E-2 

0.126  E+l 

0.1171  E-l 

-0.266  E-2 

0.226  E+l 

62.48626 

0.1507  E-2 

0.827  E-l 

0  1468  E-l 

-0.503  E-2' 

0.850 

62.99797 

0.1492  E-2 

0.166  E+l 

0.1139  E-l 

-0.334  E-2 

0.218  E+l 

15’ 

63.56852 

0.1079  E-2 

0.212  E+l 

0.1108  E-l 

-0.417  E-2 

0.196  E  +  l 

17' 

64.12778 

0.7281  E-3 

0.262  E+l 

0.1078  E-l 

-0.448  E-2 

0.200  E+l 

19' 

64.67886 

0.4601  E-3 

0.319  E+l 

0.105  E-l 

-0.515  E-2 

0.184  E+l 

21* 

65.22412 

0.2727  E-3 

0.381  E+l 

0.102  E-l 

-0.507  E-2 

0.192  E+l 

23* 

65.76474 

0.152  E-3 

0.448  E+l 

0.100  E-l 

-0.567  E-2 

0.178  E+l 

25* 

66.30195 

0.794  E-4 

0.522  E+l 

0.97  E-2 

-0.549  E-2 

0.184  E+l 

27* 

66.83663 

0.391  E-4 

0.600  E+l 

0.94  E-2 

-0.588  E-2 

0.174  E+l 

29* 

67.36933 

0.181  E-4 

0.684  E+l 

0.92  E-2 

-0.560  E-2 

0.177  E+l 

31* 

67.90051 

0.795  E-5 

0.774  E+l 

0.89  E-2 

-0.580  E-2 

0.173  E  +  l 

33* 

68.43054 

0.328  E-5 

0.869  E  +  l 

0.87  E-2 

-0.570  E-2 

0.165  E+l 

35* 

68.95972 

0.128  E-5 

0.969  E  +  l 

0.86  E-2 

-0.530  E-2 

0.174  E+l 

37* 

118.75034 

0.9341  E-3 

0.000 

0.1592  E-l 

-0.441  E-3 

0.890 

1 

that  possible  sources  of  the  excess  absorption  are 
hydrogen-bonded  dimers  and/or  clusters  of  selec¬ 
tive  size  distributions  ( 10-30  H20  molecules)  [Car¬ 
ton  and  Harden ,  1980] .  Clustering  is  assumed  to 
take  place  during  the  evaporation  process.  Another 
line  of  thought  holds  the  //-dependent  H20  at¬ 
tachment  to  hygroscopic  foreign  matter  (aerosol 


water  uptake)  responsible — a  kind  of  ‘invisible’ 
cloud  [Nilsson,  1979], 

The  nonresonant  air  spectrum  N"d  makes,  at  sea 
level  pressures,  a  small  contribution.  The  formula¬ 
tion 

N"d  =  6.2fp$3  X  10  4 
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TABLE  46. 

Data  base  for  O,  (a. 

- a ,)  and  H20  (*, ■ 

-b, )  spectral  lines  in  air  (o4 

=  a5  =  bt  =  b,  =  0). 

Center 
frequency.  v0 
(GHz) 

Strength,  a, 
(kHz/kPa) 

Temperature 
exponent,  a2 

Width,  a, 
(GHz/kPa) 

Quantum  number 
identification  (02 ) 

Lower 

Upper 

368.498350 

0.679  E-4 

0.200  E-l 

0.156  E-l 

1.  I 

2,  3 

424.763120 

0.638  E-3 

0.122  E-l 

0.147  E-l 

2,  1 

2.  3 

487.249371 

0.235  E-3 

0.122  E-l 

0.147  E-l 

2.  1 

3,  3 

715.393150 

0.996  E-4 

0.891  E-l 

0.144  E-l 

3,  3 

4.  S 

773.839732 

0.571  E-3 

0.798  E-l 

0.140  E-l 

4,  3 

4,  5 

834.145790 

0.180  E-3 

0.798  E-l 

0.140  E-l 

4,  3 

5.  5 

Identification  (H20) 

Vo 

A, 

A, 

Lower 

Upper 

•  22.235080 

0.105 

0.214  E+l 

0.281  E-l 

5,  2,  3 

6.  1,  6 

68.052 

0.180  E-2 

0.875  E+l 

0.280  E-l 

3,  2,  1  (1)“ 

4,  1,  3 

•  183.310091 

0.238  E+l 

0.653 

0.282  E-l 

2,  2,  0 

3.  1,  3 

321.225644 

0.460  E-l 

0.616  E+l 

0.220  E- 1 

9.  3,  6 

10.  2,  9 

•  325.152919 

0.155  E+l 

0.152  E+l 

0.290  E-l 

4,  2.  2 

5.  1,  5 

•  380.197372 

0.123  E+2 

0.102  E+l 

0.285  E-l 

3,  2,  1 

4,  1.  4 

386.778 

0.400  E-2 

0.733  E+l 

0.160  E-l 

11.  2,  10 

10.  3.  7 

437.34667 

0.630  E-l 

0.502  E+l 

0.150  E-l 

6.  6,  0 

7.  5.  3 

439.150812 

0.921 

0.356  E+l 

0.175  E-l 

5.  5,  0 

6,  4,  3 

443.  018295 

0.191 

0.502  E+l 

0.148  E-l 

6,  6,  1 

7.  5.  2 

•  448.001075 

0  .107  E+2 

0.137  E+l 

0.246  E-l 

3,  3.  0 

4.  2.  3 

470.888947 

0.328 

0.357  E+l 

0.181  E-l 

5.  5,  1 

6.  4,  2 

474.689127 

0.124  E+l 

0.234  E+l 

0.210  E-l 

4,  4,  0 

5.  3.  3 

488.491133 

0.256 

0.281  E  +  l 

0.222  E-l 

7.  1,  7 

6.  2,  4 

504.219 

0.380  E-l 

0.669  E+l 

0.127  E-l 

7.  7.  0 

8,  6,  3 

505.126 

0.120  E-l 

0.669  E+ 1 

0.130  E-l 

7.  7.  1 

8,  6.  2 

•  556.936002 

0.526  E+3 

0.114 

0.317  E-l 

1.0,  1 

1,  1.  0 

•  620.700807 

0.521  E+l 

0.234  E  +  l 

0.216  E-l 

4.  4.  1 

5.  3.  2 

658.340 

0.460 

0.776  E+l 

0.328  E-l 

1.0.  1(1) 

I.  1.  0 

•  752.033227 

0.259  E+3 

0.336 

0.302  E-l 

2.  0.  2 

2.  1,  1 

836.836 

0.120  E-l 

0.811  E+l 

0.170  E-l 

11.2,  9 

10.  5.  6 

859.810 

0.150  E-l 

0.799  E+l 

0.270  E-l 

2.0.  2(1) 

2.  1.  1 

899.380 

0.910  E-l 

0.784  E+l 

0.300  E-  1 

I.  1.  1  (1) 

2,  0,  2 

903.280 

0.640  E-l 

0.835  E+l 

0.280  E-l 

2.2.  1(1) 

3,  1.  2 

907.773 

0.179 

0.504  E  +  l 

0.204  E-l 

8.  3.  5 

9.  2.  8 

•  916.169 

0.890  E+l 

0.137  E  +  l 

0  249  E-l 

3.  3.  1 

4,  2.  2 

970.320 

0.940  E+l 

0.184  E+l 

0.246  E-l 

4.  3.  1 

5.  2.  4 

•  987.940 

0.145  E+3 

0.180 

0.299  E-l 

1.  1,  1 

2.  0.  2 

•  1097.368 

J.840  E+3 

0.656 

0.335  E-l 

3.0.  3 

3,  1.  2 

Read  0.307  E-3  as  0.307  x  10 
"Nonresonani  O,  spectrum  (equation  (21)). 

“Dl  denotes  doublet. 

'Rosenkranz' s  (1975]  first-order  solution  for  the  60-GHz  band  shape  has  a  shortcoming:  the  listed  values  of  a,(K  =  3  , 
5  )  have  to  be  reduced  5%  (which  is  of  negligible  consequence)  to  assure  that  a(02)  a  0.000  for  /  >  160  GHz  even  when 
N"d  =  0  (equation  (21). 

“Here  (I)  denotes  first  vibrationally  excited  slate. 


<02)  (N,) 

•  {froA/2  +  1o)J  +  2.1/td0’ x  10  7 }  ppm 

(21) 

is  taken  from  Rosenkranz  [197S]  (nonresonant  02) 
and  Stankevich  [1974]  (pressure-induced  N2).  The 


width  70  is  calculated  with  (17)  and  a,  =  0.012 
GHz/kPa  [Lam,  1977] ,  albeit  the  value  a,  =  0.0056 
generated  the  a4  and  a,  coefficients  (Table  4)  using 
Rosenkranz' s  [1975]  scheme  to  fit  laboratory  data 
[Liebe  et  al.,  1977] .  This  discrepancy  has  been 
resolved  in  a  more  profound  theoretical  treatment 
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TABLE  5.  Cloud  (fog.  haze)  attenuation  a.  in  decibels  per 
kilometer  for  1  g/m5  liquid  water  content  at  0°  and  25°C. 


Frequency  /,  GHz 

i 

10 

30  100 

200 

300 

400 

600 

800 

1000 

Equation 

(22a) 

at 

0°C  0.001 

0.09? 

0.82  5.4 

9.3 

10.7 

aw 

at 

25°C  0.001 

0.051 

0.45  4.2 

10.8 

15.3 

Equation 

(22b) 

at 

0°C 

13 

18 

23 

29 

at 

25°C 

21 

31 

40 

48 

of  the  60-GHz  band  shape  by  Smith  [1981] . 

Third,  the  liquid  water  extinction  N".  of  droplets 
with  radii  smaller  than  20  pm  (cloud,  haze,  aerosol 
hydrometeors)  can  be  derived  from  published  di¬ 
electric  data  of  bulk  water  (e\  e").  The  Rayleigh 
absorption  approximation 

N'l  =  4A9wt" /  U*'  +2)2  +  (b")2)  ppm  (22 a) 

of  Mie  scattering  losses  is  appropriate  up  to  1000 
GHz  [Falcone  et  al.,  1979] .  Frequency  and  temper¬ 
ature  dependences  of  e ',  e"  are  calculated  with  the 
Debye  model  given  by  Chang  and  fVilheit  [1979], 
which  was  found  to  describe  available  experimental 
data  up  to  300  GHz. 

At  frequencies  /  >  300  GHz,  the  rough  approx¬ 
imation 

N'i  «  0.55w/  °'a  *  ppm  (226) 

was  obtained  [Liebe,  1980]  by  fitting  e  data  reported 
for  the  submillimeter  wave  range  [e.g.,  Simpson 
et  al..  1979],  Some  numerical  examples  of  the 
attenuation  oK  =  0. 182/Af"  (in  decibels  per  kilome¬ 
ter)  are  listed  in  Table  S.  Equations  (22a)  and  (22 b) 
are  applied  to  model  the  attenuation  by  aerosols, 
haze,  clouds,  and  fog  if  their  average  water  con¬ 
centration  and  temperature  within  a  radio  path  are 
known. 

3  THE  TROPOSPHERE  PROGRAM  PI 

The  complex  refractivity  N  (equation  (2)),  as 
formulated  in  the  previous  section,  is  applied  in 
a  computer  routine  (Pi)  to  generate  values  of 
attenuation  a(/),  dispersion  £>(/),  and  refractivity 
Af„  (equation  (10)).  The  program  is  valid  for  fre¬ 
quencies  up  to  1000  GHz  and  for  heights  up  to 
30  km.  In  the  troposphere  (h  <  10  km),  only  lines 
of  H20  and  02  are  important.  Above  h  =  10  km, 
the  weak  spectral  signatures  of  trace  gases  (O,, 
CO,  N20,  etc.)  become  detectable,  but  are  not 


included  in  the  described  model.  The  calculation 
of  N(f,  P,  T,  RH)  is  done  in  a  subroutine  using 
line-by-line  superposition  at  each  frequency.  Rela¬ 
tive  humidity  RH  is  specified  and  controls,  through 

(5) ,  the  water  vapor  pressure.  A  liquid  water  con¬ 
centration  w  can  be  added.  Attenuation  a  requires 
(3),  (5),  (6),  (11),  (13),  (17),  (20a),  (21),  and  (22); 
similarly,  for  refractive  dispersion  D,  (2),  (3),  (5), 

(6) ,  (12),  (14),  and  (17)  are  applied. 


*  >,  *  '  »on  :  inr 


Fig.  I .  Attenuation  a  in  decibels  per  kilometer  and  refractive 
dispersion  D  in  parts  per  million  for  dry  (RH  =  0%)  to  saturated 
(RH  =  100%)  air  at  sea  level  (h  =  0  km)  (see  Table  3)  over 


the  frequency  range  v  =  1-1000  GHz.  Some  check  values  are 


listed  below. 

/.  GHz 

RH  =  0 

a.  dB/km 

10 

100% 

10 

0.016 

0.017 

0.028 

50 

0.299 

0.322 

0.527 

100 

0.054 

0.136 

0.886 

200 

0.028 

0.546 

5.38 

300 

0.040 

0.982 

9.71 

500 

0.108 

11.1 

116 

1000 

0.284 

124 
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Fig.  2.  Attenuation  a  in  decibels  per  kilometer  and  refractive  dispersion  D  in  parts  per  million  of  dry 
(RH  =  5%).  humid  (RH  =  50%),  and  saturated  (RH  =  100%)  air  at  sea  level  (6  =  0  km)  for  temperatures 
T  =  -5°C  (15°C— see  Figure  1)  and  35°C  over  the  frequency  range  v  =  15-350  GHz.  Continuum  model 
‘a‘  uses  equation  (20a)  and  -b'  uses  (206). 

The  main  program  PI  generates  a  frequency  grid  An  overview  to  1000  GHz  at  sea  level  is  given 
with  a  resolution  that  is  tied  to  the  line  center  in  Figure  1.  Water  vapor  attenuation  above  125 

frequencies  which  fall  within  the  desired  range.  Two  GHz  extends  up  to  3  x  104  dB/km.  The  strongest 

consecutive  values  of  v0  are  taken  from  Table  4,  maxima  of  the  rotational  H2 O  spectrum,  however, 

and  the  number  of  points  in  between  them  is  are  two  orders  of  magnitude  higher  and  fall  in  the 

specified  (usually  2-5).  Also  automatically  added  2- to  12-THz  range  [Hill  el  al.,  1980) .  Transmission 

to  the  grid  are  the  half-power  points  at  v0  ±  y.  below  25  THz  over  any  distance  is.  for  practical 

Thus  it  is  assured  that  attenuation  maxima  and  purposes,  limited  to  the  window  ranges  WI-W4 

dispersion  peaks  are  not  missed.  The  program  feeds  below  350  GHz.  The  low-frequency  wing  of  the 

one  frequency  at  a  time  into  the  subroutine  and  dispersion  spectrum,  originating  from  the  strong 

receives  values  of  a  and  D,  which  are  stored  557-GHzH201ine,  makes  its  influence  known  down 

temporarily  in  two  files.  When  all  numbers  of  the  to  about  100  GHz. 

grid  have  been  processed,  the  output  is  presented  Atmospheric  EHF  (30-300  GHz)  transfer  proper- 
in  a  numerical  printout  and  processed  by  a  graphics  ties  are  exemplified  in  Figure  2  at  sea  level  for 

routine.  Examples  are  depicted  in  Figures  1-4  taking  the  temperatures  -5°  and  +35°C  and  the  difference 

*v  =  0.  Some  obvious  conclusions  drawn  from  these  between  the  two  water  vapor  continuum  spectra 

figures  are  as  follows.  (20a)  and  (206)  is  illustrated.  The  transparency  of 
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Fig.  3.  Attenuation  a  in  decibels  per  kilometer  and  refractive  dispersion  D  in  parts  per  million  Tor  the  oxygen 
microwave  band  at  the  altitudes  h  =  0-30  km  (see  Table  3)  covering  the  frequency  range  i>  =  50-70  GHz. 


the  window  ranges  W2-W4  is  dominated  by  the 
controversial  water  vapor  excess  absorption,  which 
is  normally  modeled  by  (20a).  The  temperature 
dependence  of  window  attenuation,  expressed  in 
terms  of  e ,  is  only  T  ',  but  changes  to  T'1  when 
RH  is  the  governing  variable.  A  more  detailed 
picture  of  molecular  attenuation  appears  in  Figure 
3.  The  oxygen  microwave  spectrum  dominates  in 
the  range  50-70  GHz.  Close  to  sea  level  the  60-GHz 
lines  are  merged  into  an  unstructured  band  shape, 
the  maximum  intensity  of  which  is  pressure-propor¬ 
tional  until  the  lines  separate  (h  >  15  km).  Above 
IS  km,  the  overlap  disappears  and  radio  channels 
with  up  to  400-MHz  bandwidth  can  be  accommo¬ 
dated  between  the  lines.  Details  of  the  center  portion 
of  the  60-GHz  band  are  magnified  in  Figure  4.  This 
figure  shows  interesting  dispersion  properties, 
which  change  with  increasing  height  from  negative 
to  positive  gradients  dD/dv  between  the  5*  and 
5  lines.  Numerical  values  depicted  in  Figures  1, 
2,  and  4  may  serve  as  control  numbers  for  Program 
PI. 

4.  THE  MESOSPHERE  PROGRAM  P2 
Program  P2  calculates  the  attenuation  of  single 
lines  as  they  appear  above  h  =  30  km.  The  formula¬ 


tion  starts  with  an  isolated,  pressure-broadened  line 
described  by  (3)  and  (11)  for  N"c  =  0,  (15),  (16), 
and  (17).  Decreasing  pressure  in  (17)  eventually 
approaches  the  finite  Doppler  line  width 

y„  =  6.20 kHz  (2.3) 

where  m  is  the  molecular  weight  (e.g.,  02,  m  = 
32,  0  =  I,  v0  =  60  GHz;  yD  =  65.8  kHz).  As 
a  consequence,  a  change  to  a  Gaussian  line  shape 
function  takes  place. 

The  convolution  of  Lorentzian  and  Gaussian 
shape  functions  is  called  the  Voigt  profile,  which 
is  governed  by  the  parameter  y  =  y/yD  and  which 
is  appropriate  when  this  ratio  falls  in  the  range 
®  1 0-0. 1 .  Numerical  evaluation  of  the  height-depen- 
dent,  complex  Voigt  function  requires  considerable 
computational  effort.  Usually,  the  Voigt  profile 
‘frequencies’  x  and  ‘pressures’  y  are  normalized 
in  multiples  of  y0  and  the  intensities  are  expressed 
with  reference  to  the  Lorentzian  shape.  Frequency 
profiles  reduced  in  such  a  manner  are  labeled  for 
absorption  w(x,  y)  =  F" /  F\  and  for  dispersion  v(x, 
y)  =  F' IF",  .  Examples  are  shown  in  Figure  5 
together  with  the  pressure  profiles  of  the  two 
maxima  u0  and  v0  approaching  zero  intensity.  Maxi- 


ATTENUATION  AND  PHASE  OF  RADIO  WAVES 


1193 


I 


.197 

r . 

■  .153 


8.14 


/5+ 


»  *T  *  ’ 


9.57 

4.99 


-J 


- 

3.43 

/  *\ 

/2.91H8 

/  f  j]  \ 

- 

2.84  \ 

/  ' ;  '  \ 

i  7+ 

z 

^  3 

«r/  , 

Z 

.>6 

1 

1  V  1 

- 

1  -57 

*T-"r-  T’ 

60.0 

f'RCOUENCi  -GHZ 


.126 

r 

.099^ 

Kn 


T>  -?97 


T// 

i  ,  / 


!! 1 

U  i 

vv127 

-.132  •u' 


/ 

\  : 

’  Curve 

h 

RH 

2 

J 

km 

X 

4  ^ 

:  1 

0 

50 

.032 

:  2 

10 

5 

■* 

:  3 

20 

1 

4 

30 

0 

S9.6  60.:.’  60.1 

I'REOUOC)  -  'GHZ 


V,-, 


Fig.  4.  Attenuation  a  in  decibels  per  kilometer  and  refractive  dispersion  D  in  parts  per  million  at  the 
oxygen  microwave  band  center  for  the  altitudes  h  =  0-30  km  (see  Table  3)  covering  the  frequency  range 
v  =  59.4-60.7  GHz. 


mum  line  attenuation  is  followed  by  P2  through 
the  mesophere  until  it  has  dropped  below  0.01 
dB/km.  A  FORTRAN  IV  algorithm  was  adopted 
for  the  calculations  of  m(jc,  y)  as  reported  by 
Pierluissi  et  al.  [1977] .  An  adequate  approximation 
can  be  made  by  replacing  -y  in  (17)  with 

•Y*  =  V%J  +Y»  (24) 

and  keeping  the  shape  functions  (equations  (15)). 
Isolated  line  attenuation  is  then  simply  ((3),  (11), 
(15).  (16),  and  (24))  given  by 

6, (/)  =  “.,/(!+  z2)  dB/km  (25) 

where  a„  =  0.1820 fS/yh  is  the  maximum  at  i>„. 
and  z  =  (u,,  -  f)/y„  is  the  normalized  frequency. 

Zeeman  splitting  of  isolated  oxygen  lines  due  to 
the  influence  of  the  earth’s  magnetic  Held  strength 
H  introduces  considerable  complications  [Lenoir. 
1968] .  The  presence  of  the  steady  field  H  (in  units 
of  gauss)  splits  each  of  the  K  ’  lines  into  three 
groups  of  (2K  ’  ±  1)  sublines,  thereby  redistributing 


Fig.  5.  Normalized  Voigt  profiles  | Faddryea  and  Tarentev. 
1961)  for  absorption  u(x,  v)  and  dispersion  vix.  y)  Shown  are 
frequency  profiles  (v  =  /)  at  p  =  const  (solid  curves)  and 
pressure  profiles  of  maximum  absorption  u„y  and  peak  disper¬ 
sion  |  (dashed  curves).  For  Oj  lines,  y"  =  a,8"\  and  for 
H;0  lines,  y"  =  A,  ft"". 
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Fig.  6.  Schematic  energy  level  diagram  displaying  the  Zeeman 
components  for  K  =  3  ‘  and  3  oxygen  microwave  lines.  K 
is  total  rotational  angular  momentum  quantum  number  in  odd 
integers.  J  -  K  and  K  ±  I  is  the  rotational  quantum  number. 
M  is  the  magnetic  quantum  number  (Table  6),  and  g  is  the 
Lande  splitting  factor. 


the  line  attenuation  a,  (equation  (25))  over  a  fixed 
frequency  range. 

Part  of  program  P2  is  a  subroutine  which  calcu¬ 
lates  the  line  center  frequencies 

V',  =  V„  +  •>)(*:.  M)28.03  x  10  4H  GHz  (26) 

and  the  line  strength  coefficients 

a]  =  a,£.(A.  M)  kHz/kPa  (27) 

for  the  Zeeman  components  of  each  K  '  microwave 
line.  The  components  are  identified  by  the  magnetic 
quantum  number  M.  The  rotational  quantum 
number  J,  which  is  determined  by  the  total  quantum 
number  K,  sets  the  limits  on  Af.  Figure  6  gives 
an  example  of  the  schematic  distribution  of  energy 
levels  and  allowed  transitions  for  the  case  K  = 
3.  The  Lande  g  factor  is  the  constant  of  propor¬ 
tionality  for  a  splitting  of  each  J  level.  The  fact 
that  g  is  different  for  levels  J  >  1  leads  to  the 
anomalous  Zeeman  effect,  while  the  two  lines  K 
=  1 '  display  the  normal  Zeeman  effect.  The  cal¬ 
culation  procedure  for  the  relative  frequency  shifts 
i)(/C  Af)  and  relative  intensities  §( AT,  Af)  is  sum¬ 
marized  in  Table  6.  The  formulations  are  based 
on  Lenoir's  [1968]  work  with  two  exceptions:  (1) 
the  intensities  £(K,  Af )  are  normalized  to  the  theo¬ 
retical  definition  of  the  coefficient  a ,  [Liebe  and 
Gimmestad.  1978] .  and  (2)  both  K  lines  are  referred, 
against  convention  for  K  ,  to  the  same  set  of  Af 
values  (see  Figure  6),  which  leads  to  the  limits  of 
Af  specified  in  Table  6.  An  02  line  breaks  up  into 
numerous  components,  which  are  organized  in  three 
groups  labeled  it.  o  * ,  and  cr  . 


TABLE  6  Relative  frequency  shift  q  ( M .  K)  and  relative  intensity  factor  i(M.  K  )  for  Zeeman  components  of  O,  microwave 

lines 


Attenuation  Zeeman 


K  "  line 


rate 


transitions  q(Af.  K )  <  ±  I 
«<l  -  K ) 


UM.  K)  <  1 
3((A  +  IE  -  M' | 


U(M.  K )  <  +1 
M{K  +  2) 


K  line 

UM.  K)  -  I 
3( Af  4  l)(A  ;  -  M: ) 


(AM  =  0) 


(AM  =  I) 


(AM  =  -I) 


K(K  4  I)  (K  4  l)(2Af  4  1  )(2  A'  4  3)  A  (g  4  I)  A(2A  4  I  )<2A :  4  K  -  1 1 

M  limits:  +  M  «  A  A  -  I  >  M  =-  -(A  -  I) 

M(l  -  A)  -  A  3(A  4  M  4  IHA  4  M  4  2)  (M  4  IRA  4  2)  -  I  3(A  4  IRA  -  M«A  -  M  -  I) 


A(A4l)  4(A  4  l)(2A  4  IK2A  4  3)  A(A4|)  4A(2A  4  l«2AJ  4  A  -  I) 

M  limits:  ±M«A  A-2>M?-A 

M(I-A)4A  3<A  -  M  4  IXA  -  M  4  2)  (M  -  IHA  4  2)  4  I  3(A  4  IMA  4  MM  A  4  M  -  I) 


A(A4l)  4(A  4  l)(2A  4  l)(2A  4  3)  A(A  4  I)  4A(2A  4  l)(2A"  4  A  -  I) 

M  limits:  ±MsA  A»M»-(A-2) 


...A 
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Frequency  Deviation  MHz 


Fig.  7.  Zeeman  attenuation  patterns  of  the  oxygen  microwave  line  K  =  3'  for  altitudes  h  =  30-100  km  (see  Table  3).  Each 
frame  displays  ir(«,).  (t'(u,).  and  o  («,)  patterns  for  the  magnetic  field  strengths  //  =  0.3  (left  hand)  and  H  =  0.6  G  (right 
hand).  The  patterns  are  symmetric  with  respect  to  the  center  axis  (interchange  <?'  and  tr  ).  The  frequency  deviation  An  = 
/  -  v„  is  between  -*40  MHz  for  h  =  30  and  ±2  MHz  for  h  =  60-100  km.  The  maximum  attenuation  rale  q„  is  that  of 
the  isolated  unsplit  line  (//  =  0);  the  value  in  parentheses  is  calculated  with  program  PI.  The  values  of  in  decibels 

per  kilometer  are  for  0.3  (see  Figure  9)  and  0.6  G 


In  principle,  three  different  Zeeman  patterns  al2, 
are  possible  for  any  K  ’  line.  These  patterns  are 
obtained  by  evaluating  all  of  the  17  and  4  coefficients 
in  the  tr  and  a 1  groups  and  applying  a  line  shape 
such  as  (25).  Each  sum  of  sublines  determines  a 
Zeeman  pattern  according  to 

« I.M  =  k(M,  K)/{l  +  [(/-  v„(M,  K))/yh\ !) 

MK 

(28) 

where 

a„  =  O.I82f(alp/yJd’  exp  (a2(l  -  8)]  dB/km 

The  unsplit  line  H  =  0  follows  from  (28)  by  assuming 
4(A#,  K)  =  1  and  v2  =  i>„.  For  more  accurate 


calculations,  the  Lorentzian  (equation  (15))  in  (28) 
may  be  replaced  by  the  Voigt  profile  (Figure  5). 

Additional  information  has  to  be  provided  before 
the  mesospheric  02  line  pattern  a .  can  be  calculated, 
that  is,  the  polarization  of  the  radio  wave  and  its 
orientation  with  respect  to  the  geomagnetic  field 
H.  For  linearly  polarized  radiation  an  angle  of 
orientation  0  is  defined  between  the  wave's  magnetic 
field  component  in  a  plane  of  constant  phase  (i.e.. 
perpendicular  to  the  direction  of  propagation)  and 
the  geomagnetic  field  direction  H;  for  circularly 
polarized  radiation  the  angle  6  is  defined  between 
the  plane  of  constant  phase  and  H.  A  mesospheric 
03  line  a,  then  generally  consists  of  a  mixture  of 
the  three  Zeeman  patterns: 


Fig  X  Same  as  Figure  7.  except  K  -  19  oxygen  microwave  line. 

Wave  polarization  Zeeman  pattern  a  . 

linear  a,  cos’  6  +  (a,  +  a,)  sin’ t) 

circular  right-handed  2a,  sin-’  9  +  0.5(a,  +  a,  +  a,)  cos'  9  (29) 

circulat.  left-handed  2a,  sin:  9  +  0.5(a,  +  a_,  +  a,)  cos:  9 


Equations  (29)  imply  that  mesospheric  02  line 
attenuation  is  polarization  and  direction  (anisotropic 
medium)  dependent.  Further  coordinate  trans¬ 
formation  is  necessary  when  a  fixed  antenna  re¬ 
ceives  radiation  from  a  mesospheric  space  element. 

The  patterns  aI2,  of  the  02  lines  K  =  1 ±  to 
29 e  have  been  calculated  for  local p-T-H  conditions 
over  the  height  range  30-100  km  at  two  magnetic 
field  strengths,  H  =  0.3  and  0.6  G.  Two  examples, 
K  =  3  *  and  19  ,  are  depicted  in  Figures  7  and 
8.  (The  complete  set  K  1  =  1-29  is  available  upon 
request  (to  be  published  in  a  report).)  Above  h 
=  70  km  the  individual,  now  mostly  Doppler- 
broadened  components  become  discernable.  Al¬ 
though  the  details  depend  on  specific  K  values. 


the  overall  features  at  given  heights  are  similar  to 
the  cases  illustrated  in  Figures  7  and  8  except  for 
different  scale  factors. 

The  maximum  of  the  pattern  a ,  ( v„ )  for  h  =  30-100 
km  is  shown  in  Figure  9  covering  K  =  1 ±  to  29 ±, 
so  far  as  a  threshold  of  0.02  dB/km  is  exceeded. 
The  K*  =  19  line  pair  has  a  values  independent 
of  height  (temperature)  over  the  range  h  =  20-50 
km.  The  results,  in  combination  with  a  beginning 
dissociation  to  atomic  oxygen  (e.g.,  at  h  -  120 
km,  02/0  =  0.25),  establish  h  =  100  km  as  the 
plausible  boundary  to  outer  space  for  radio  path 
modeling.  The  cutoff  height  for  water  vapor  is  much 
lower  ( h  =  20  km). 

Having  a  full  Zeeman  picture  at  hand,  one  might 
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ZtCHII  Effect 


Fig.  9.  Maximum  attenuation  a,(v„,  H  =  0.3  G)  for  the  oxygen  microwave  lines.  K  -  I  *  to  29 1  (Table 
4a)  at  altitudes  h  =  30-100  km  (see  Table  3). 


consider  applications  for  this  information.  The  uni¬ 
form  mixing  of  02  affords  opportunities  to  remotely 
sense  mesospheric  temperature  profiles  T(h)  or 
geomagnetic  field  strength  including  variations,  H 
±  dH,  via  emmissions  originating  from  a .(h).  A 
multiplet  pattern  is  expanded,  with  increasing  H, 
like  an  accordion  (see  Figure  7),  at  the  most  to 
vn  ±  2.5  MHz  (equation  (26)).  Variations  in  the 
geomagnetic  Field  strength  H  translate  into  height 
dependent  changes  of  selective  attenuation  a.[h). 
The  example  in  Figure  10  reveals  that  the  response 
h(dH)  for  0|(v(|)  is  most  pronounced  in  the  region 
h  =  60-70  km. 

Power  saturation  of  a  mesospheric  02  Zeeman 


pattern  is  our  final  concern.  An  active  system  (i.e., 
radar),  in  principle,  can  diminish  or  even  eliminate 
the  calculated  attenuation  a..  A  rough  estimate  for 
the  onset  of  saturation  effects  yields  that  the  radiat¬ 
ed  power  density  at  v0  needs  to  exceed  values  on 
the  order  of  0.5p2  in  units  of  W  cm'3  Pa  2,  which 
is  not  out  of  the  question  above  h  =  70  km. 

5.  CONCLUSION 

Moist  air  was  characterized  as  an  atmospheric 
propagation  medium.  The  basic  physics  of  molecu¬ 
lar  absorption  in  a  radio  path  have  been  cast  into 
a  model  with  optimum  computer  run  time,  but 
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Fig.  10.  Response  8  of  the  maximum  attenuation  a,(v„)  in 
percent  to  an  equal  percentage  drop  in  the  magnetic  field  strength 
H  over  the  altitude  range  h  =  40-100  km  (see  Table  3)  for 
the  oxygen  microwave  lines  K  =  1  ’  to  29  ‘ . 


without  undue  approximations.  The  paper  has 
expanded  upon  an  EHF  version  (30-300  GHz)  by 
Liebe  and  Gimmestad  [  1978]  in  terms  of  frequency 
and  height  coverages  and  has  provided  graphical 
and  numerical  examples.  In  addition,  the  formula¬ 
tion  of  dispersion  has  been  clarified.  Shortcomings 
of  the  present  model  lie  in  the  omitted  spectroscopic 
data  base  for  trace  gases,  and  in  the  lack  of  a 
' — model  that  describes  the  /(//-dependent  water  up**'' 
takeTf*-^  w)  by  various  atmospheric  aerosol  dis¬ 
tributions  [Nilsson,  1979] . 

Use  of  programs  PI  and  P2  is  straightforward. 
The  routines  (available  in  FORTRAN  IV  upon  re¬ 
quest)  are  operated,  even  by  nonexpert  users, 
simply  by  calling  the  program  name,  choosing  from 
a  battery  of  model  atmospheres,  and  typing  in  the 
interactively  requested  input  parameters  within  the 
limits  of  (9):  meteorological  data  are  converted  into 
propagation  data.  Detrimental  effects  such  as  fluc¬ 


tuations  in  amplitude,  phase,  and  group  velocities, 
and  the  direction  of  the  radio  wave  may  be  predicted 
from  spatial  and  temporal  variations  of  P-RH-T. 
A  priori  assessments  can  be  made  of  rare  events 
for  reliable  system  operations  and  of  possible  band¬ 
width  limitations. 

The  model  may  serve  as  a  reference  for  compari¬ 
son  with  current  research.  This  is  especially  valuable 
with  respect  to  studies  of  water  vapor  excess 
absorption,  which  is  accounted  for  only  by  an 
empirical  correction.  A  reliable  expression  for  re- 
fractivity  N(f,  P,  RH,  T,  H,  6)  would  be  a  boon 
to  workers  in  the  field  since  N  (equation  (2))  is 
at  the  hub  of  all  millimeter  and  submillimeter 
propagation  problems  through  the  nonprecipitating 
atmosphere.  An  important  task  ahead  is  research 
into  the  nature  of  H20  excess  absorption,  already 
underway  in  several  laboratories  including  ours  at 
the  Institute  for  Telecommunication  Sciences, 
NTIA. 
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